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Pair correlations between large transverse momentum neutral pion triggers (pr = 4-7 GeV/c) and 
charged hadron partners (j>t = 3-7 GeV/c) in central (0-20%) and midcentral (20-60%) Au+Au 
collisions are presented as a function of trigger orientation with respect to the reaction plane. The 
particles are at larger momentum than where jet shape modifications have been observed, and the 
correlations are sensitive to the energy loss of partons traveling through hot dense matter. An out- 
of-plane trigger particle produces only 26 ± 20% of the away-side pairs that are observed opposite of 
an in-plane trigger particle. In contrast, near-side jet fragments are consistent with no suppression 
or dependence on trigger orientation with respect to the reaction plane. These observations are 
qualitatively consistent with a picture of little near-side parton energy loss either due to surface bias 
or fluctuations and increased away-side parton energy loss due to a long path through the medium. 
The away-side suppression as a function of reaction-plane angle is shown to be sensitive to both the 
energy loss mechanism in and the space-time evolution of heavy-ion collisions. 

PACS numbers: 25.75.Nq,25.75.Bh 



I. INTRODUCTION 

Collisions of heavy nuclei at the Relativistic Heavy 
Ion Collider have created matter with energy densities 
exceeding the predicted threshold for deconfinement of 
color charge into a hot dense plasma In this quark 
gluon plasma(QGP), quarks and gluons are not bound 
within hadronic states and the matter behaves collec- 
tively. Comparisons with hydrodynamic simulations in- 
dicate rapid thermalization of the colliding system into 
a hot dense nuclear medium. The produced medium af- 
fords an opportunity to study the properties of a new 
phase of quantum chromodynamics (QCD) in an extreme 
environment. 

Hard scattering with large momentum exchange be- 
tween partons in the incoming nuclei is well-described 
by perturbative QCD (pQCD). The scattered partons 
emerge back-to-back in azimuth in the plane transverse 
to the beam direction, and fragment into a pair of cor- 
related cones of high momentum particles, referred to as 
jets. The study of jets and their hadronic fragments in 
heavy-ion collisions provides insight into the properties of 
hot dense nuclear matter. Measurements of single high 
transverse momentum (pr) particles Q and correlations 
between high-p T particles 0t3 have demonstrated that 
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the fast partons embedded in the produced medium dis- 
sipate a large amount of their initial kinetic energy. 

In this paper, we present angular correlations of hadron 
pairs with both hadrons in the midrapidity range < 
0.35. Fragments from the same jet form a peak at small 
relative azimuthal angle (A</>), i.e. the near-side peak. 
Pairs composed of one fragment from each jet will ap- 
pear in an away-side peak at A(f> ~ ir. Past measure- 
ments 0-01 for hadrons > 5 GeV/c have shown that the 
away-side correlations peak is suppressed relative to base- 
line measurements in p+p collisions. The suppression of 
the away-side jet is a signature of parton energy loss in- 
side the medium. The same measurements show that 
near-side jet fragments at large momentum are not sup- 
pressed. This feature of the data is understood to result 
from the requirement of a large momentum particle in 
the final state, which creates a bias towards small energy 
loss either by the preferential selection of hard scatter- 
ings near the medium surface [6| or due to fluctuations 
in energy loss Q. 

The detailed mechanism by which partons lose en- 
ergy when passing through a deconfmed medium are 
not yet fully understood. In pQCD descriptions of the 
parton-medium interaction the predicted parton energy 
loss should scale as the path length squared [8|]. In 
competing anti-de-Sitter space/conformal field theory de- 
scriptions characterizing a strongly coupled medium, the 
energy loss scales as the path length cubed |8| . The vari- 
ation in azimuthal angle of the away-side jet suppression 
with respect to the reaction plane ((j) s ) is sensitive to 
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the total amount of energy lost by the away-side parton 
along long paths (out-of-plane) or short paths (in-plane) 
through the medium. The degree to which the away-side 
jet suppression varies will be determined in part by the 
path-length scale of energy loss. 

Single particle observables at high py, such as the 
nuclear suppression with respect to the reaction plane 
(Raa(4>s)) or the azimuthal anisotropy (i.e. 1)%) |9j, are 
also sensitive to this path-length variation of energy loss. 
Current pQCD calculations predict a lower V2 than is 
found in the data and may imply a larger than path 
length squared dependence to energy loss [H| • The 
reaction-plane dependence of the back-to-back jets pro- 
vides an additional test on the path length dependence 
in that the two particle observable selects a different dis- 
tribution of hard-scattering locations and should probe 
longer paths through the medium than single particle 
observables. The path-length dependence of both single 
and two particle observables have already been studied 
through selection of the collision centrality [J, Q . How- 
ever, centrality selection varies not only the path length, 
but also other important properties of the medium (e.g. 
the overall energy density). Selection with respect to 
the reaction plane more directly varies the path lengths, 
while leaving the other medium properties unchanged. 

In addition to the uncertainties associated with the en- 
ergy loss mechanisms, many of the details within hydro- 
dynamic simulations of heavy- ion collisions have not been 
fully constrained and tested by experiment. For instance, 
one such uncertainty is the geometrical description of the 
energy deposited by the colliding nuclei which could con- 
tribute to the degree of away-side suppression variation 
with respect to the reaction plane. Two competing de- 
scriptions, the Glauber model [l2[ and the Color Glass 
Condensate [l3| , predict different azimuthal distributions 
of matter with respect to the reaction plane. Thus the 
two descriptions give different in-plane and out-of-plane 
path lengths through the medium. These descriptions 
are also used as different starting points to the hydro- 
dynamic evolution of the medium. Other model uncer- 
tainties include, but are not limited to, the extent of ge- 
ometry fluctuations, the time required for thermalization 
into a hydrodynamic medium, the characteristics of the 
phase transition to confined hadrons, and the conditions 
under which those hadrons become free-streaming parti- 
cles into the vacuum. These ambiguities in the proper 
modeling of heavy-ion collisions can result in significant 
uncertainty in the extracted properties of the medium, 
such as the shear viscosity 

In midcentral collisions (the middle 20-60% of the total 
cross section) the variation of the away-side suppression 
is expected to be largest as the collision zone is the most 
anisotropic. In contrast, central collisions are much more 
isotropic and so provide a sample of events with small 
anisotropy expected in the away-side suppression. For 
instance, in the Glauber model, midcentral events will 
have a root-mean-square thickness through the medium 
of 3.2 fm in the in-plane direction versus 4.8 fm in the 



out-of-plane direction, which is a 50% variation in path 
length. However, for central 0-20% collisions, the path 
length through the medium varies from 5.0 fm in the in- 
plane direction to 5.8 fm in the out-of-plane direction, 
which is a much smaller 16% variation. It is notable that 
the thickness through the medium in midcentral colli- 
sions changes more with respect to the reaction plane 
than it does between central and midcentral collisions 
where the away-side suppression at large momentum is 
already known to vary [J] . Also worth noting is that the 
largest thickness in midcentral collisions is comparable to 
the shortest thickness in central collisions. 

Any prediction for the away-side suppression with re- 
spect to the reaction plane will be a convolution of the 
energy loss and a description of the space-time evolu- 
tion of the medium. In the limit where the medium is 
never fully opaque to fast partons and the energy lost by 
the typical parton is some fraction of its initial energy, 
the away-side suppression will increase with angle with 
respect to the reaction plane. This results because the 
average path length through the medium of the recoil 
parton is longer when out of the reaction plane. It is 
possible, in the extreme limit of a medium with a large 
opaque core and thin transparent corona, the away-side 
suppression could instead weaken as the trigger particle 
orientation varies from the in-plane to out-of-plane direc- 
tions. The weakening in the thin corona scenario results 
from two effects; a larger relative number of scattering 
centers producing a pair of back-to-back final-state par- 
ticles in the out-of-plane direction, but also the variation 
of the trigger particle multiplicity by angle with respect 
to the reaction plane. However, it is worth noting that 
a large core and thin corona is an extreme configuration. 
Variations within more realistic models of the away-side 
suppression will be intermediate between these extreme 
scenarios. 

In this paper, we present azimuthal correlation mea- 
surements between pairs of neutral pion trigger particles 
(t) within pj, = 4-7 GeV/c and charged hadron associ- 
ated partner particles (a) within p r = 3-4, 4-5, and 5-7 
GeV/c. These combinations of final-state particle mo- 
mentum ranges have previously been shown to be dom- 
inated by jet fragmentation as they are above medium- 
induced two-particle correlations which contribute sig- 
nificantly at lower momenta [3, § . The low momentum 
structures (the "ridge" and "shoulder" ) may b e the result 
of parton- medium interactions (e.g. [15l4l8| |) or global 
correlations from fluctuating initial conditions (l9l [20l j . 
These fluctuations have substantially less impact at large 
pair momentum where the background contribution be- 
comes small. In this study, as illustrated in Fig. [TJ us- 
ing only particles at large pair momentum, the away-side 
(A</> ~ tt) suppression by trigger particle orientation with 
respect to the reaction plane (<p s = (j> — is presented 
as a probe of both the mechanism of parton energy loss 
and the space-time evolution of matter created by the 
collisions of large nuclei. 

A previous measurement [2l| by the STAR coUabo- 
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FIG. 1: (Color online) Definition of azimuthal angles. Trigger 
and associated partner particles are measured at <\? and 4> a , 
respectively. The trigger particle orientations are taken with 
respect to the reaction plane, 4> s — — if). The relative 
azimuthal separation of the trigger particle and the associated 
partner particle is = (j> — (f> a . 



ration for 20-60% centrality between 4-6 GeV/c trigger 
particles and 2 < p^ < p f T GeV/c partner particles for 
two 45° wide in-plane and out-of-plane selections indi- 
cated an increased suppression of the away-side jet for 
the out-of-plane, but with little significance due to large 
underlying event subtraction uncertainties. The new re- 
sults presented in our paper have sufficient statistics to 
specify a trend in the away-side suppression in midcentral 
collisions at larger momentum where subtraction uncer- 
tainties are negligible. 



II. EXPERIMENT 

The results presented here are based on 3.4 billion 
minimum-bias Au+Au events recorded by the PHENIX 
detector in 2007. Comparisons to p+p collisions use pre- 
viously published measurements from data recorded in 
2006 [5]. Collision centrality was determined by divi- 
sion into percentile of the integrated charge collected by 
beam-beam counters (BBC) [231 located at \rj\ between 
3.0 and 3.9. The timing between the arrival of charged 
particles in the north and south BBC was used to re- 
construct the event position along the collision axis (z- 
vertex), and to restrict the event sample to ±30 cm of 
the nominal interaction point of the two beams. The 
orientation of reaction-plane azimuthal angle {ip) is re- 
constructed event-by-event using the quadrupole compo- 
nent (^2) of the charge in the Reaction Plane (RXPN) 
detector [23|, located at \q\ between 1.0 and 2.8. The res- 
olution of the RXPN detector is highest for midcentral 
collisions (~20%) where both the quadrupole component 
and the detector occupancy are large. The set of resolu- 
tion corrections, A„ : n € {2,4,6,8}, for single particle 



anisotropics, v n , where: 



are estimated from correlations between the independent 
north (i/jn) and south (ips) RXPN reaction-plane recon- 
structions [3, Hfl . A single fit parameter (a;) is mapped 
into the resolution corrections via: 




The fit parameter is extracted from the correlations via: 



+z(l (z) + M*))+y (h(z) + L 1 (z)) j(3) 
where 

x 2 

z = — cos (ip N - ip s ) (4) 

The set of functions, J2fc and L 2 k, are the even-ordered 
modified Bessel functions and the modified Struve func- 
tions respectively. The extracted values used to correct 
the measured second-order azimuthal anisotropy, are A2 
= 0.66(4) and 0.66(3) for 0-20% and 20-60% collisions, 
respectively. A 10% systematic uncertainty in 0-5% col- 
lisions and 5% elsewhere accounts for non-flow contri- 
butions to the resolution corrections (23[. The similar 
values are a result of the peak in reaction-plane resolu- 
tion appearing near 20% centrality. A direct inspection 
of these reaction-plane distributions for events contain- 
ing a photon above 1 GeV/c did not reveal significant 
contributions from jets. 

Neutral pion trigger particles are reconstructed from 
photon clusters measured by either lead-glass or lead- 
scintillator electromagnetic calorimeters (EMCal) in the 
two central arms of PHENIX, in total covering \rj\ < 0.35 
and 2 x 90° in azimuth [2y|. Clusters are subject to cuts 
based on the known response of the EMCal, including 
noisy and low-response towers, as well as shower shape 
cuts. Neutral pions are identified through the 27 de- 
cay channel by pairing all photons within an event. In- 
correct pairings between photons create a broad com- 
binatorial background under the ir° mass peak. This 
background is minimized by requiring the reconstructed 
mass to lie near the 7r° mass peak. This requirement was 
0.125-0.160 MeV/c 2 for central events, but was relaxed 
to 0.120-0.165 MeV/c 2 in midcentral events where the 
combinatorial background is lower. Since combinatorial 
pairs are more often made with the abundant photons 
found at low energy, the energy asymmetry of the decay 
(\Ei - E 2 \/{E 1 + E 2 )) was restricted to be less than 0.5 
for 0-5% central events. This was also relaxed slowly for 
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more peripheral events until all pairs with asymmetries 
less than 0.7 were accepted. The tightness of the cuts 
was used to control the rate of combinatorial pairings 
such that 7T° trigger particles have a signal-to-background 
ratio averaged over the mass window of 4:1 in central col- 
lisions and 10:1 in midcentral collisions. 

Charged hadron partner particles are reconstructed in 
the central arms using the drift chambers (DC) with hit 
association requirements in two layers of multi-wire pro- 
portional chambers with pad readout (PCI and PC3), 
achieving a momentum resolution, Ap/p, of 0.7%©l.l%p 
(GeV/c). Only tracks with unambiguous and distin- 
guishable DC and PCI hit information are used. Pro- 
jections of these tracks are required to match a PC3 hit 
within a ±2a proximity window to reduce background 
from conversion and decay products. A track association 
to a signal in the Ring Imaging Cerenkov detector is used 
to reject electrons for partner selections below 5 GeV/c 
where little signal is produced by charged pions. 



III. PAIR ANALYSIS 

Within an event, all pairs formed from ir° trigger par- 
ticles (pj, = 4-7 GeV/c) and three sets of charged hadron 
associated partner particles {pj, = 3-4, 4-5, 5-7 GeV/c) 
are measured. Two centrality classes are used: a central 
selection of 0-20% collisions and a midcentral selection 
of 20-60% collisions. Trigger particles are separated into 
six 15° bins in azimuthal angle with respect to the re- 
action plane, <f) s — (f> — ij). The angular resolution of 
the measured reaction plane, at approximately 25°, is 
larger than this binning; consequently, significant smear- 
ing takes place between neighboring trigger orientation 
bins. Pairs within PHENIX are collected at different ef- 
ficiencies due to the non-uniform central arm acceptance. 
The relative pair efficiencies are corrected by mixed pair 
distributions in which trigger and partner particles are 
drawn from different events of the same class (bins of 
5% centrality, 5 cm z-vertex). The resulting acceptance- 
corrected distributions are reported as correlation func- 
tions, C(A(j)), which are defined as: 



dm'" 



C(A0) 



dAtj, 



dAtj, 



■dAcj) 



dAtj, 



I 



dmt 



dAtj, 



■dA<j> 



(5) 



where m' a is the number of measured pairs per event for 
either the same or mixed events, as indicated. Double- 
struck notation (n) is used here to indicate measured 
quantities. Representative correlation functions for in- 
plane and out-of-plane trigger particle orientations are 
shown in Fig. [5] The full set of the measured correla- 
tion functions used in this analysis is shown in Figs. 1131 
and[T5] Note that these distributions are not corrected 
for reaction-plane resolution. 

These inclusive pairs are assumed to correlate in one 
of two ways. (1) Two particles within the same event 
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FIG. 2: (Color online) Correlation functions for the most in- 
plane, <^ s =0-15°, (solid squares) and out-of-plane, <^ 3 =75-90°, 
(open squares) trigger 7r° orientations in central 0-20% and 
midcentral 20-60% collisions, left and right columns respec- 
tively, for 3-4 GeV/c partner hadrons. Expected underlying 
event contributions are shown as solid curves (see text for 
details). 



may correlate trivially by participation in the same col- 
lision geometry. These pairs produce an azimuthal an- 
gular correlation from the the single particle anisotropy 
with respect to the reaction plane. (2) Two particles 
may also correlate with each other via the same hard- 
scattering process. These particles will be fragments from 
the same (di)jet. To separate the jet particle pairs from 
the other back grou nd pairs, the two-source assumption 
is expressed as [27j : 



C(A<f>) = J(A0) 

+ b ( 1 + -cos(2A0) + -cos(4A0) ) (6) 

\ a a J 

where the jet contribution to the correlation function is 
contained in J(A(f>). The remaining harmonic terms de- 
scribe the background contribution which is complicated 
by the trigger particle binning with respect to the re- 
action plane. The background modulation coefficients 
(a, /3, 7) are calculated via: 



„ t . . . sin (2c) . 
l + 2v* cos(2<ft fl ) o v ' A 2 
2c 

2v\ cos (40 s ) A 4 

2v\v a 2 + 2v a 2 (1 + v\) cos (20 s ) A 2 



(7) 



2c 



„ t n , , , \ sin (4c) . 
2vlv 2 cos (4<(> s ) — A — L A 4 
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„ „ i . ( . sin (6c) . 
+ 2«>*cos(60 s ) ; ; A 6 

DC 

7 = 2v 4 v 4 I + 2<(l + ^)cos(40 s ) 



(8) 



sin (4c) 
4c 



A, 



+ 2u|«^ cos(2^ s ) 



sin (2 c) 



2c 



A, 



. , . sin (6c) 
cos(60 s ) „ v ; A 6 



6c 



2W4W4 cos ( 



sin (8c) 



(9) 



This description of the background accounts for the trig- 
ger particle binning and reaction-plane resolution effects 
on the background shape [28|. The trigger particle ori- 
entation appears explicitly in terms of the bin center, <f> s , 
and width, c. Single particle anisotropy values, V2 and 
i>4, were measured by correlating the trigger and partner 
particles with respect to the reaction plane, such that: 

C(0-V) = 1 + 2v° bs cos(2(</> - V)) 

+ 2v° hs cos(4(</> - V)) (10) 

where the observed anisotropies are corrected for the 
reaction-plane resolution as described previously in 
Equation [1] Given sufficient detector resolution and nar- 
rowness of the trigger particle orientation binning, the 
sign of the cos(2A(/>) term in Eq. [B]will flip sign between 
in-plane and out-of-plane bins as shown in Fig.[2j This ef- 
fect is expected as selecting out-of-plane trigger particles 
should decrease the likelihood of finding a second back- 
ground particle nearby. The same is not true for parti- 
cles correlated via hard scattering. Both the second- and 
fourth-order anisotropy of the background correlations 
have been considered as the finite fourth-order contribu- 
tions were determined to be non-negligible for some trig- 
ger particle orientations. Likewise, higher-order terms in 
the reaction-plane resolution correction are also included. 

The uncertainties on the reaction-plane resolution cor- 
rections (A n ) and the observed anisotropies (i>2 bs i v 4 hs ) 
are propagated separately as they impact the away-side 
suppression with respect to the reaction plane in char- 
acteristically different ways. The uncertainty in the 
reaction-plane resolution corrections is fully correlated 
between trigger orientations. For instance, this uncer- 
tainty increases (or decreases) both the extracted in- 
plane and out-of-plane jet yields at A<j> = n. However, 
the uncertainty in the observed anisotropies is fully anti- 
correlated between trigger orientations. Thus, this un- 
certainty increases the extracted in-plane yield while de- 
creasing the out-of-plane yield (or vice versa). At large 
momentum, both of these subtraction uncertainties are 
small and always dominated by other sources. 

The subtraction procedure was also examined for con- 
tamination of the jet correlations by fakes in the charged 
tracking, which become significant at large pt- The fake 
high pt tracks are present only in the partner sample and 
are largely uncorrelated with trigger particle for the part- 
ner pt presented here. Thus the fake tracks, which are 
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FIG. 3: Correlation functions for trigger 7r°s averaged over 
all trigger orientations in central 0-20% and midcentral 20- 
60% collisions, left and right columns respectively, for 3-4 
GeV/c partner hadrons. Expected average underlying event 
contributions are shown as solid curves. 



already less influential in events with a high pt 7r°, are 
subtracted with other uncorrelated pairs as part of the 
background contribution, so long as the anisotropies are 
measured with the same particle cuts. The subtraction 
robustness against tracking fakes at high pt was checked 
by repeating the procedure with a 3a PC3 matching re- 
quirement. 

By taking the trigger particle orientation as <j) s = 7r/4, 
the bin width as c = 7r/2, and by truncating higher 
than second-order terms, the functional form of the back- 
ground in Eq. [6] reduces to the v\ x modulation used in 
previous trigger particle orientation averaged results such 
as those found in This property is demonstrated in 
Fig. [3] where the trigger particles from all orientations are 
considered. 

The background level, bo, is determined using the zero 
yield at minimum (ZYAM) method [13]. At high-p T , the 
well-separated near- and away-side jets provide a large 
angular region at mid-A</> angles with negligible jet sig- 
nal. This allows the ZYAM level to be found with neg- 
ligible bias and sufficient statistics despite the lower ef- 
ficiency PHEN1X has for collecting pairs near 90°. The 
ZYAM uncertainty was estimated through simulation of 
the statistical uncertainties as has been described in [2^] . 

The jet contribution, J(Acj)), is then reported as a per- 
trigger azimuthal yield such that: 



1 



1 



dAcj) 



J dAcj) 



J(A0). 



(11) 



The efficiency-corrected single particle and pair rates are 
ciency. 



n respectively. The single particle partner effi- 
e a , is estimated in simulations of detector accep- 
tance and occupancy as was done in By design, the 
trigger particle efficiency cancels in the ratio. 



IV. RESULTS 

Central events, 0-20% collisions, are analyzed as a 
cross check against experimental artifacts in midcentral 
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collisions since they have a smaller away-side jet yield. 
Thus, the central events should exhibit a smaller trigger 
particle angle variation, require a larger reaction-plane 
resolution correction, a larger event correlation subtrac- 
tion, and have increased background in ir° identification. 
Representative per-trigger azimuthal yields in central col- 
lisions for each of the partner momentum selections for 
the most in-plane and most out-of-plane trigger particle 
selections are shown in Fig. |4j Figure [TT] shows the full 
set of the measured per-trigger azimuthal yields used in 
this analysis for central collisions. The most in-plane and 
most out-of-plane trigger-particle orientations select the 
shortest and longest average path lengths through the 
medium, respectively, and thus may be expected to have 
the maximimum differences. 

On the near-side, a jet distribution is clearly observed 
for each selection. A direct comparison between the most 
in-plane and most out-of-plane trigger shows no signif- 
icant variation. The measurement at mid-A0 demon- 
strates the good agreement resulting from correct de- 
scription of the underlying event correlations. On the 
away-side, the jet yield is small due to medium sup- 
pression and the statistical precision suffers once the 
pairs are divided among the various trigger particle ori- 
entations. No evidence of experimental artifacts such 
as over-subtraction or incorrect description of the back- 
ground is seen, despite the challenging analysis environ- 
ment present in central collisions. 

Integrated near- and away-side per-trigger yields (Y) 
are calculated within angular Acj) windows, as indicated 
in Fig. |H approximating the 2a width of the jet distri- 
butions measured in the trigger particle orientation av- 
eraged results. The near-side azimuthal integration win- 
dows are Acj) < tt/9 (< 3tt/18) for p^ > 4 GeV/c (< 4 
GeV/c). Similarly, the away-side azimuthal integrations 
windows are n-A<f> < 3ir/l8 (< 2tt/9) for > 4 GeV/c 
(< 4 GeV/c). Use of these windows corresponds to an 
assumption that the jet distributions do not widen sig- 
nificantly at high pt, as a function of the trigger particle 
orientation with respect to the reaction plane. This as- 
sumption is supported by the absence of significant cen- 
trality dependence in jet correlation widths (< 20%) for 
particles at high pr @- Within statistical uncertainties 
a constant jet width is consistent with the data. Inte- 
grated yields as a function of trigger particle orientation 
for both the near- and away-side are then corrected for 
the reaction-plane resolution. The resolution correction 
is applied such that: 



l + 2(v° hs ' Y /A 2 ) cos(20 s ) 
1 + 2^2 bs ' Y cos (2(j) s ) 



mcas V ys 



)• (12) 



where Y and Y melLS are the corrected and uncorrected 
yields, respectively. The value of is the observed 

second-order anisotropy of integrated per-trigger yield 
with respect to the reaction plane and is determined by 
fitting the trigger particle orientation dependence of each 



,) measurement individually. This procedure is 



the similar to the correction of reaction-plane resolution 
on single particles, here applied to integrated per-trigger 
pair yields. 

The corrected per-trigger yields (Y) are reported as 
the nuclear jet suppression with respect to p+p collisions, 
Iaa = Ya+a/Y p+p . The result for central collisions is 
shown in Fig. [5] The variation of the fit used in the reso- 
lution correction is the dominant source of (/> s -correlated 
uncertainty, having larger impact than the insignificant 
event anisotropy uncertainties. In the case of zero signal 
variation with reaction plane orientation, the correction 
becomes completely correlated with statistical scatter in 
the uncorrected measurement. Thus, the ^-correlated 
systematic uncertainty from the resolution correction is 
conservatively treated as correlated with the statistical 
uncertainty when computing the final significance of the 
measured trends. For the same reason, this source of 
systematic uncertainty has little correlation between the 
centrality and momentum selections. 

For central events the near-side suppression is consis- 
tent with a constant as a function of cj> s within the sta- 
tistical and ^-correlated systematic uncertainties. The 
values are also consistent with no suppression when con- 
sidering the global scale uncertainty that appears on the 
trigger particle orientation averaged Jaa- On the away- 
side, there is significant suppression in central events, as 
evidenced by the trigger particle averaged Jaa, but the 
statistical precision with which to determine the cj> s vari- 
ation is limited. 

Mid-central events, 20-60% collisions, have greater ec- 
centricity and could be expected to show correspond- 
ingly larger trigger particle orientation dependence due 
to path-length variation through the collision zone. The 
same set of representative per-trigger azimuthal yields is 
shown in Fig. [5] for the midcentral selection. The full set 
at midcentrality is shown in Fig. 1121 Again, the near-side 
jets for the most in-plane and most out-of-plane trigger 
particle orientations are consistent with each other, a di- 
rect indication of little variation with respect to the re- 
action plane. The mid-A</> are also in agreement with 
zero, as before, further demonstrating that the under- 
lying event flow correlations are well described by Equa- 
tions[B][ni In contrast to the near-side, the away-side mea- 
surements (see insets in Fig. [6]) change between the in- 
plane and out-of-plane trigger particle orientations with 
the latter having consistently smaller yield for all partner 
momenta. 

The integrated near- and away-side per-trigger jet 
yields for midcentral collisions are shown in Fig. [Jj The 
near-side jet is essentially flat, with negligible suppres- 
sion (Jaa (0s) = !)• The away-side jet yield is increas- 
ingly suppressed with increasing cf> s . This falling trend 
results in only small associated particle yield remaining 
for out-of-plane trigger particle orientations. 

In order to quantify the variation and significance of 
the trigger particle orientation dependencies shown in 
Figs. [5] and \7\ the ratio of the out-of-plane to in-plane 
suppression (^aa/^aa) i s constructed. In the ratio, the 
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FIG. 4: (Color online) Per-trigger azimuthal jet yields for 
the most in-plane, S =O-15° (solid circles) and out-of-plane, 
<f>s =75-90° (open circles) trigger particle selections in central 
0-20% collisions for various partner momenta. Insets show 
away-side region on a zoomed scale. Bars indicate statistical 
uncertainties. Underlying event modulation systematic un- 
certainties are represented by bands through the points while 
the corresponding normalization uncertainties are shown as 
dashed lines around zero. Near- and away-side jet yield inte- 
gration windows are indicated with arrows. 
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FIG. 5: (Color online) Nuclear jet suppression factor, Jaa, 
by angle with respect to the reaction plane, <f>„, for near- and 
away-side angular selections, circles and squares respectively, 
in central 0-20% collisions for various partner momenta. Bars 
indicate statistical uncertainties. The shaded band shows the 
systematic uncertainty on the reaction-plane resolution un- 
smearing correction. Solid points show trigger particle angle 
averaged results and the global scale uncertainty. 



global scale uncertainties on each measurement cancel. 
The Taa values at S = 0° (4" A ) and at 90° are es- 

timated by both linear and flow-like cosine fits to the trig- 
ger particle angle measurements and evaluation at these 
angles. The reported ratios are therefore independent of 
the chosen binning with respect to the reaction plane and 
the values do not rely heavily on the assumed functional 
form of the dependence. The best-fit was determined by 



X 2 minimization in which: 

~2 _ (Vi j[ e sys (J sys,i ~ / (0s)) ,2 /i o\ 

X — ~2 ( f \ "+" 6 sys K 1 ^) 

where e sys is ±1 for the ±la sys variation of the (j) s - 
correlated systematic error [3(| ■ As discussed above, the 
systematic uncertainty is conservatively treated as fully 
correlated with the statistical uncertainty. The differ- 
ence between linear and cosine fits provides only a small 
source of additional uncertainty due to the unknown true 
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FIG. 6: (Color online) Per-trigger azimuthal jet yields for 
the most in-plane, <^ 3 =0-15°, (solid circles) and out-of-plane, 
<f>s =75-90° , (open circles) trigger particle selections in mid- 
central 20-60% collisions for various partner momenta. In- 
sets show away-side region on a zoomed scale. Bars indicate 
statistical uncertainties. Underlying event modulation sys- 
tematic uncertainties are represented by bands through the 
points while the corresponding normalization uncertainties 
are shown as dashed lines around zero. Near- and away-side 
jet yield integration windows are indicated with arrows. 
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FIG. 7: (Color online) Nuclear jet suppression factor, Jaa, 
by angle with respect to the reaction plane, <f> 3 , for near- 
and away-side angular selections, circles and squares respec- 
tively, in midcentral 20-60% collisions for various partner mo- 
menta. Bars indicate statistical uncertainties. The shaded 
band shows the systematic uncertainty on the reaction-plane 
resolution unsmearing correction. Solid points show trigger 
particle angle averaged results and the global scale uncer- 
tainty. 



functional form. The resulting values of ^aa/^aa an d 
the total uncertainty are shown in Figs. [8] and [9] The 
average value of Iaa/ITa across partner momentum is 
constructed by weighting the individual measurements 
by the p+p per-trigger yields @ . In general, the data are 
well fit by both the linear and cosine functions, giving rea- 
sonable x 2 - No evidence is seen for systematic deviations 
from either fit within the sizable statistical uncertainties 
and both forms give similar goodness of fit values. These 
values appear along with the / A a / -^aa ratios in Table Q] 
For both central and midcentral collisions, the near- 



side jet yield is independent of trigger particle orientation 
with respect to the reaction plane within one standard 
deviation of the experimental uncertainties. These mea- 
surements are consistent with surface bias of the hard 
scattering center created by the requirement of a trigger 
particle and a resulting short path length through the 
collision zone traversed by the near-side parton. Cen- 
tral collisions have insufficient statistics to determine the 
away-side variation. 

However in midcentral collisions where the expecta- 
tion of surface bias would lead to a large variation in 
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TABLE I: Angle with respect to the reaction-plane dependence of the nuclear suppression factor, 7aa, expressed as the ratio 
between in-plane and out-of-plane trigger particles from linear and cosine fits to the data (see text for details). The total 
uncertainty taking into account the correlations between sources is reported. 
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FIG. 8: (Color online) Angle with respect to the reaction- 
plane dependence of the nuclear suppression factor, Iaa, ex- 
pressed as the ratio between in-plane and out-of-plane trigger 
particles from fits to the data in central 0-20% collisions. The 
bars represent total uncertainty taking into account the cor- 
relations between sources (see text for details). 



FIG. 9: (Color online) Angle with respect to the reaction- 
plane dependence of the nuclear suppression factor, Iaa, ex- 
pressed as the ratio between in-plane and out-of-plane trigger 
particles from fits to the data in midcentral 20-60% collisions. 
The bars represent total uncertainty taking into account the 
correlations between sources (see text for details). 



the path length traversed by the away-side parton, the 
measurements show a significant falling trend with in- 
creasing trigger particle angle with respect to the reac- 
tion plane. The suppression of away-side jet fragments in 
the out-of-plane direction is larger than in the in-plane 
direction, the out-of-plane away-side jet peak having only 
(26 ±20)% of the yield of the in-plane direction. Thus the 
large variation by angle with respect to the reaction plane 
is significant. Assuming the modulation to be flow-like 
(dominated by the second-order variation), the suppres- 
sion pattern implies «2 AA — 0.29j^*f . As the midcentral 
away-side measurements are consistent between p^ selec- 
tions within the stated uncertainties, the hint of a rising 
trend in p9p is not significant. The values quoted here are 
consistent with those previously measured in [2~lj and 
provide a factor four better constraint in the Jaa/^aa 
ratio. 

Recent single particle measurements of azimuthal 



anisotropy at high px (6 — 9 GeV/c) found that V2 = 
0.13 ± 0.01 ± 0.01 If. Thus, the away-side per-trigger 
yields at high pt favor an anisotropy larger than that 
measured for the single particles. However, the difference 
is marginal and additional measurements will be needed 
to confirm. 

Shown in Fig. [10] are the results of a Monte-Carlo en- 
ergy loss calculation from T. Renk [3l],[32| using the time- 
space evolution provided by two different hydrodynamic 
simulations (33l . l34j and two initial state descriptions, 
Glauber and CGC. These particular combinations of a 
jet energy loss model and collision evolution together pre- 
dict less variation in the away-side suppression with re- 
spect to the reaction plane than is witnessed by the data. 
Variation of the initial geometry description within [33[ 
between Glauber and CGC produces only small changes 
in the extracted Iaa out-of-plane to in-plane ratio, in- 
dicating limited sensitivity to this model parameter of 
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FIG. 10: (Color online) Away-side J A a74a ratio for mid- 
central, 20-60% collisions, from Fig. [5] The solid line shows 
the results from an energy loss calculation [3ll . l32l ] using two 
hydrodynamic evolution models 33, 34]. The shaded band 
shows the uncertainty that results from the selection of a par- 
ticular hydrodynamic evolution; the lower extent covering [33j 
and the upper covering [34|. Dotted lines show the uncer- 
tainty from the initial event geometry (Glauber or CGC) as 
calculated within 1331. 



the reaction plane dependent dijet observable. However, 
other model parameters that vary between the two hy- 
drodynamic models (such as the thermalization time and 
freeze-out temperature) were found to impact the away- 
side suppression anisotropy to a greater degree, indicat- 
ing sensitivity to simulation parameters that are not well- 
constrained by other measurements. Consequently, these 
data warrant more detailed study with various energy 
loss models, and also different space-time evolution mod- 
els. 



V. SUMMARY 

We have shown that away-side jet fragment suppres- 
sion increases substantially with increasing angle with 
respect to the reaction plane in midcentral Au+Au col- 
lisions at ys^" = 200 GeV. The away-side yield in the 
out-of-plane orientation is reduced by a factor of ~ 4 rel- 
ative to the in-plane direction. In contrast, the measured 
near-side Iaa is reaction plane independent, and consis- 
tent with no suppression. These results directly show 
that the energy lost by fast partons in the hot nuclear 
medium increases as their paths through the medium 
become long. A theoretical description of these exper- 
imental data implementing an energy loss formalism and 
a time-space evolution of the collision should be sought in 
union with other experiment al q uantities; such as Raa, 
7 A a, and i? AA (0 s ) H S S SI [H. Energy loss for- 
malisms that have successfully described the large mo- 
mentum i?AA and ^aa may be paired with a particular 
time-space evolution in also describing the <f) s dependence 
of these same quantities. As shown for the combination 



above, the data presented here disagree with the present 
calculations. These data should play an important role 
in constraining simulations of the space-time evolution 
of heavy-ion collisions and the subsequent extraction of 
medium properties. 
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Appendix: Appendix 

The complete set of per-trigger yields and correlation 
functions used as source material for the analysis of the 
dependence of the away-side suppression on angle with 
respect to the reaction plane are shown in Figs. [TT1[T11 
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FIG. 11: (Color online) Per-trigger azimuthal jet yields in central 0-20% collision for all trigger orientations with respect to 
the reaction plane in 15° selections from the most in-plane to the most out-of-plane proceeding from left to right. Associated 
partner momentum selections from lower to higher p^ are arranged from top to bottom. Bars indicate statistical uncertainties. 
Flow modulation systematic uncertainties are represented by bands through the points while flow normalization uncertainties 
are shown as dashed lines around zero. 
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FIG. 12: (Color online) Per-trigger azimuthal jet yields in midcentral 20-60% collision for all trigger orientations with respect 
to the reaction plane in 15° selections from the most in-plane to the most out-of-plane proceeding from left to right. Associated 
partner momentum selections from lower to higher p^ are arranged from top to bottom. Bars indicate statistical uncertainties. 
Flow modulation systematic uncertainties are represented by bands through the points while flow normalization uncertainties 
are shown as dashed lines around zero. 
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FIG. 13: (Color online) Correlation functions in 0-20% central collisions for all trigger orientations with respect to the reaction 
plane in 15° selections from most in-plane to most out-of-plane proceeding from left to right. Associated partner momentum 
selections from lower to higher arranged from top to bottom. Expected flow contributions are shown as solid curves (see 
text for details). 
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FIG. 14: (Color online) Correlation functions in 20-60% central collisions for all trigger orientations with respect to the reaction 
plane in 15° selections from most in-plane to most out-of-plane proceeding from left to right. Associated partner momentum 
selections from lower to higher arranged from top to bottom. Expected flow contributions are shown as solid curves (see 
text for details). 



15 



[1] K. Adcox et al., Nucl. Phys. A757, 184 (2005). [17 

[2] S. S. Adler et al., Phys. Rev. C76, 034904 (2007). [18 
[3] J. Adams et al., Phys. Rev. Lett. 97, 162301 (2006). 

[4] A. Adare et al., Phys. Rev. C78, 014901 (2008). [19 

[5] A. Adare et al., Phys. Rev. Lett. 104, 252301 (2010). [20 

[6] C. Loizides, Eur. Phys. J. C49, 339 (2007). [21 

[7] T. Renk and K. J. Eskola, Phys. Rev. C 75, 054910 [22 

(2007). [23 

[8] F. Dominguez, C. Marquet, A. H. Mueller, B. Wu, and [24 

B.-W. Xiao, Nucl. Phys. A811, 197 (2008). [25 

[9] A. Adare et al. (2010), 1006.3740. [26 

[10] S. A. Bass, C. Gale, A. Majumder, C. Nonaka, G.-Y. [27 

Qin, T. Renk, and J. Ruppert, Phys. Rev. C 79, 024901 [28 

(2009). 

[11] S. Wicks, W. Horowitz, M. Djordjevic, and M. Gyulassy, [29 

Nucl. Phys. A784, 426 (2007). 

[12] M. L. Miller, K. Reygers, S. J. Sanders, and P. Steinberg, [30 

Ann. Rev. Nucl. Part. Sci. 57, 205 (2007). [31 

[13] H.-J. Drescher and Y. Nara, Phys. Rev. C 75, 034905 [32 

(2007) . [33 
[14] M. Luzum and P. Romatschke, Phys. Rev. C78, 034915 

(2008) . [34 
[15] J. Casalderrey-Solana, E. V. Shuryak, and D. Teaney 

(2006), hep-ph/0602183. [35 
[16] S. S. Gubser, S. S. Pufu, and A. Yarom, Phys. Rev. Lett. 
100, 012301 (2008). 



I. Vitev, Phys. Lett. B630, 78 (2005). 

A. D. Polosa and C. A. Salgado, Phys. Rev. C75, 041901 

(2007). 

P. Sorensen, arXiv:0808.0503 [nucl-ex]. 
J. Takahashi et al. (2009), 0902.4870. 
J. Adams et al., Phys. Rev. Lett. 93, 252301 (2004). 
M. Allen et al., Nucl. Inst, and Meth., A 499, 549 (2003). 
A. Adare et al., Phys. Rev. Lett. 105, 062301 (2010). 
S. Afanasiev et al., Phys. Rev. C80, 024909 (2009). 
S. Afanasiev et al., Phys. Rev. C80, 054907 (2009). 
A. Adare et al., Phys. Rev. Lett. 101, 232301 (2008). 
S. S. Adler et al., Phys. Rev. Lett. 97, 052301 (2006). 
J. Bielcikova, S. Esumi, K. Filimonov, S. Voloshin, and 
J. P. Wurm, Phys. Rev. C 69, 021901 (2004). 
A. Sickles, M. P. McCumber, and A. Adare, Phys. Rev. 
C81, 014908 (2010). 

A. Adare et al., Phys. Rev. C77, 024912 (2008). 
T. Renk, Phys. Rev. C78, 034904 (2008). 
Renk, T., private communication. 

C. Nonaka and S. A. Bass, Phys. Rev. C 75, 014902 
(2007). 

K. J. Eskola, H. Honkanen, H. Niemi, P. V. Ruuskanen, 
and S. S. Rasanen, Phys. Rev. C 72, 044904 (2005). 
A. Adare et al., Phys. Rev. C80, 024908 (2009). 



